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Abstract 

We consider an extension of the (t—U) Hubbard model taking into 
account new interactions between the numbers of up and down elec- 
trons. We confine ourselves to a one-dimensional open chain with L 
sites (4 L states) and derive the effective Hamiltonian in the strong 
repulsion (large U) regime. This Hamiltonian acts on 3 L states. We 
show that the spectrum of the latter Hamiltonian (not the degenera- 
cies) coincides with the spectrum of the anisotropic Heisenberg chain 
(XXZ model) in the presence of a Z field (2 L states). The wave func- 
tions of the 3^-state system are obtained explicitly from those of the 
2 L -state system, and the degeneracies can be understood in terms of 
irreducible representations of U q (sl(2)). 
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1 Introduction 



The one-dimensional Hubbard model, which depends on two parameters (U 
and t), has received a lot of attention due to its possible application (in its 
2-d version) to describe high-T c oxide materials [1] . The large U limit of this 
system was discussed in Refs. [2]. 

In this paper we consider the following Hamiltonian: 



H = H + Hi , 
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Here c+ CT , c, ; 



and ni iCr = c^ ct q i0 -((t =t, I) are creation, annihilation and oc- 
cupation number fermionic operators. In the Hubbard model one has the 
hopping term (with coefficient —1), the repulsion term (with coefficient ^) 
and the external field (1.2b) (coefficient ^). We have generalized this Hamil- 
tonian by adding a quartic interaction term (coefficient j) and an external 
field type term (coefficient |£). The term proportional to ^ corresponds to 
a boundary term; it does not change the thermodynamical properties of the 
system but, as we shall see, plays an important role in its symmetry proper- 
ties. For similar reasons we have added the constant y, which will be shown 
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to be useful for the identification of representations of some Hecke algebras. 
We have no phenomenological reasons to introduce the supplementary terms 
in the Hamiltonian but, as will be seen, the phase structure of the physical 
system thus obtained is much richer. 

If, as in the t-J model [1], we take the large j limit keeping ^ = g, j = 
v i ft = w i Yt = a an< ^ T = ^ fi xe d, as shown in Appendix A, we obtain instead 
of a Hamiltonian acting on 4 L states, an effective Hamiltonian acting on 3 L 
states since the remaining states get an infinite energy and drop out of the 
spectrum. This effective Hamiltonian reads 

H' = H' Q + H[, (1.3) 

where 

K = £ Ui (i.4) 

i=l 

and 

Ui= - (g~ ot+i + qTq7+i + T t T i+ 1 + T i T t i) 
+ + «;(e? + 6^) 

+ atf-elJ + b (1.5) 



and 



L 



^ = -»E e *- ( L6 ) 



i=i 



We use here the notation of Ref . [3] : 
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(1.7) 

Voooy Voioy 

Obviously: 

[Hi S z ] = [H' , S + ] = [H' , S-] = [H'„ S°] = , (1.8) 

where 

s* = £E^^ = E e f;s° = I>?. (L9) 

Z i=l i=l i=l 

Thus if£ is 5C/(2) <g> C/(l) invariant (generators 5 ± ,5 2 and 5°). As will 
be seen later, the symmetry algebra is much larger, namely H' Q is invariant 
under U q (sl(2)) <g> U(l) for any value of q ! 

Special cases of the quantum chain (1.4) were already obtained in a 
mathematical context. In Ref. [4] it was noticed that if 

v — , w + a = q , w — a = q' 1 ,6 = (1-10) 

the C/j's are generators of the Hecke algebra. We notice that if we take: 

v = 0,w = l, a = (1.11) 

we obtain a special case of the one-parameter-dependent U q (sl(2))(q 3 = 1) 
invariant quantum chain of Ref. [3] . In this model one uses nilpotent repre- 
sentations. The same quantum chain plays an important role in a physical 
application [5]. If one considers the one-dimensional problem in which two 
types of molecules, say A and B, diffuse in a lattice (vacancies represent the 
third state of the problem) and undergo chemical reactions, the time evolu- 
tion operator that appears in the master equation is given, for certain rates, 
by the Hamiltonian (1.4). We stress that in order to obtain Eq. (1.4) the 
molecules A and B have to play a symmetric role. 

It is the aim of this paper to study the spectrum and eigenf unctions of 
the Hamiltonian H' Q . If this problem is solved, since H[ commutes with H' Q , 
the spectrum of H' is trivially obtained for finite chains. The calculation in 
the thermodynamical limit of the spectrum and correlation functions in the 
presence of H[ has to be done, as usual, with care, since the ground state of 
the system is dependent on the coupling constant g appearing in Eq. (1.6). 
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As will be shown in the present paper, the spectrum and wave functions 
of the Hamiltonian H' Q can be obtained from those of another Hamiltonian, 
namely the XXZ spin \ Heisenberg chain in the presence of a Z field, which 
is integrable. The last Hamiltonian depends (leaving aside the boundary- 
terms) on the anisotropy A and the strength h of the Z field. These two 
parameters are related to the two parameters v and w in Eq. (1.5). It turns 
out that the energy levels of the Heisenberg chain (2 L states) and those of 
the chain given by Eq. (1.4) (3 L states) are identical. The degeneracies of 
the 3 L states chain depend on those of the 2 L chain, which in turn depend 
on the values of A, h and boundary terms. 

The paper is organized as follows. In Sec. 2 we shortly review the phase 
diagram of the Heisenberg chain since the Hamiltonian (1.4) has the same 
phase diagram. We also remind the reader of the known symmetries of the 
model (quantum group symmetries) and point to an amusing new symmetry, 
which has as consequence the degeneracy of just two special energy levels 
for any finite length (see also Appendix B). We then rewrite the Heisenberg 
model in a different basis and show how to get other quantum chains with 
the same spectrum, the Hamiltonian (1.4) being the simplest case. 

In Sec. 3 we give a simple construction that allows, starting from an 
eigenfunction of the Heisenberg chain corresponding to a given eigenvalue, 
to get all the corresponding eigenfunctions of the 3 L states chain with the 
same energy. 

In Sec. 4 we look at this problem from a different point of view. Starting 
from the observation that H' is invariant under the SU (2) given by the gen- 
erators and S z , we try to g-deform the chain and look for a Hamiltonian 
that commutes with the known generators and S z of the quantum algebra 
U q (sl(2)) (for q = 1, coincides with S^). Surprisingly enough, the chain 
stays g-independent although the generators are (using the coproduct 
rules) g-dependent. This implies that the chain is invariant under U q (sl(2)) 
for any q, in particular it is invariant under the affine U q (sl(2)) algebra for 
any q. In U q (sl(2)) one combines U q (sl(2)) with U q -i(sl(2)). This allows us to 
understand the degeneracies of the (1.4) chain in terms of finite-dimensional 
irreducible representations of U q (sl(2)). 

Starting from this observation, in Sec. 5 (see also Appendix C) we ask 
under which circumstances a 3 L chain can be invariant under U q (sl(2)), if 
we take a spin and a spin | on each site. The answer we get is again the 
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Hamiltonian (1.4). We warn the reader that, as opposed to the first sections, 
the last two (4 and 5) are rather mathematical in character. 

Finally, in Sec. 6 we describe the degeneracies of the 3 L chain if the 
underlying Heisenberg chain has supplementary symmetries as described in 
Sec. 2 and Appendix B. Our conclusions are presented in Sec. 7. 

2 A procedure to obtain quantum chains with 
the same spectrum as the Heisenberg chain 

The spin | anisotropic Heisenberg chain is defined by the Hamiltonian 



where 

v, = -i 

2 



afaf +1 + o\a\ +x + A(T?<rf +1 + h(a z + a z i+1 ) + a{a\ - <j? +1 ) + (3 

(2-2) 

and a x , o y and o z are Pauli matrices. As is well known, this chain is integrable 
with the help of the Bethe ansatz. If A, h and a are arbitrary, there are no 
degeneracies in the spectrum. Obviously, 

[I',H"\=0, (2.3) 

where 

I Z = \jl°l- (2-4) 
z i=i 

One obtains higher symmetries in three cases: 
Case la 

A = 0, h + a = -q, h-a = -g" 1 , (3 = -(q - q' 1 ) . (2.5) 

In this case the Hamiltonian is C/ g (sZ(l/l))-symmetric [7] and its diagonalized 
form can be expressed [8] through fermionic occupation numbers: 

H = e n + J2 (q + " 2 cos (^j ^ n k , (2.6) 
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where eo = is a zero fermionic mode and the n^'s (i = 0, • • • , L — 1) take 
the values and 1. Barring accidental degeneracies, each energy level of the 
Hamiltonian (2.6) is twice degenerate due to the zero fermionic mode as long 
as q is real. If q is a root of unity, additional degeneracies occur (see Ref. 
[8]). The V^'s defined by Eqs. (2.2), (2.5) are the generators of the Hecke 
algebra: 

VMttVi -Vi = v t±1 v^ l±1 - v i±1 

[Vi, V}] = , >2 (2.7) 

V^iq + q-'m (i = l,...,L-l) 
corresponding to the quotient [9]: 

(y i V l+2 )V l+1 {q + q- l -V i ){q + q- 1 -V l+2 ) = Q (i = 1, 2, . . . , L - 3) . (2.8) 

Case lb 

q + q~ l q-q' 1 n q + q' 1 , „s 

A = — - — ,h = 0,a= — - — ,/3= — . (2.9) 

The Hamiltonian is f/ g (s/(2))-symmetric [10], the V^'s are again the genera- 
tors of the Hecke algebra (2.7) but now correspond to the Temperley-Lieb 
quotient: 

ViVi^Vi = Vi . (2.10) 

Cases 2 and 3 

(A±h) 2 = l + a 2 . (2.11) 

This symmetry is a bit unusual. It just says that the unique energy level 
with l z — \ has the same energy as one level with I z — |f — 1. Alternatively 
the level with I z = — ^ has the same energy as one level with I z — — ^ + 1. 
For details see Appendix B. 

At this point we will make a change of basis in the Hamiltonian (2.1) 
suggested by the quantum chain formulation of diffusion-reaction processes 
done in Ref. [5]. For an iV-state model (N = 2 for the Heisenberg chain, iV = 
3 for the Hamiltonian (1.4)) we will denote the states by m — 0, 1, . . . , AT — 1. 
The state m = will play a special role. In the space of iV x N matrices 
we take the usual basis of matrices E mn (m, n — 0, 1, . . . , N — 1). The matrix 
E mn has as only no n- vanishing element the one in the m row and n 
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column, and its value is 1. In this basis the Pauli matrix cr z can be written 

as 

<j z = E 00 — E u = 1 — 2E 11 . (2.12) 
In the new basis the V^'s of Eq. (2.2) have the following expression: 

+ w (El 1 + EH,) + « W - E lli) + b > (2-13) 

where 

v = -2A, w = A + h, a = a, b = - + hj • (2.14) 

In order to construct an iV-state system with the same spectrum, we 
write non-diagonal terms that do not mix the states m — 1,2, . . . , N — 1 and 
diagonal terms where one permutes the states "m" (the state "0" excluded). 
This takes us to the following Hamiltonian: 

H=J2V t , (2.15) 

i=i 



where 



JV-l 

Vi = - £ {E^E^ + E?°E°™) + vey i+l 

m=l 

+ w (e° + e° i+1 ) + a (e° - e° +1 ) + b (2.16) 



and 



JV-l 
\ ^ rpmm 



E mm . (2.17) 

m=l 

In particular, if iV = 3, we recover the Uj's of Eq. (1.5). Since the states 

"0" play a special role it is useful to introduce an operator Z, which counts 
them: 

Z = X>°° = L -X>° , (2-18) 

i=i i=i 

its eigenvalues z are 0, 1, . . . , L. In particular if N = 2: 

Z = ^ + /\ (2.19) 
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i.e., it coincides with the magnetization. 

In Sec. 3 we will show, for N >2, that the spectrum of the Hamiltonian 
(2.15) is independent of the number of states N. One corollary of this ob- 
servation (which has also been checked by direct calculation) is that the V^'s 
of Eq. (2.16), if the conditions (2.5) or (2.9) are satisfied, also satisfy the 
Hecke algebra (2.7) and the corresponding quotient relations (2.8) or (2.10). 
We have to stress that all these results are valid for open chains only. The 
chains with periodic boundary conditions have spectra which depend on the 
number of states N (see Sec. 3). 

At this point it is useful to remind the reader of the phase diagram [11] 
of the Heisenberg chain, which is shown in Fig. 1. For h = 0, the system is 
massive with a ferromagnetic ground state if A > 1, massless and conformally 
invariant if —1 < A < 1, and again massive with an antiferromagnetic 
ground state if A < — 1. For a given h, the system is massive ferromagnetic 
if A > 1 — h, then undergoes a Pokrovsky-Talapov [12] transition at A = 
Apt = 1 — h, is in a massless incommensurate phase for A c < A < 1 — h, 
and reaches again the massive antiferromagnetic phase at A < A c . The 
A c = A c (h) curve is obtained by solving the equation 



Notice that the A = A PT line corresponds to the special symmetry related to 
the condition (2.11) for a — 0. The same phase diagram is valid through the 
translation given by Eq. (2.14) for all the iV-state systems described by Eqs. 
(2.15) and (2.16), and therefore for the infinite repulsion limit of the extended 
Hubbard model. In this case, in order to interpret the phase diagram, it is 
better to work with a fixed value of the electron nearest-neighbour interaction 
(v) and to vary by changing the fugacity (w). Phase (III) corresponds to a 
phase where the fermions and holes are physically separated. In phase (I) 
there are no fermions in the ground state, and in phase (II) we have fermions 
and holes mixed. 




(2.20) 



9 



3 Energy levels and wave functions of the N- 
state model 



In this section, we write the wave functions of the Hamiltonian (2.15), (2.16) 
in terms of those of the Hamiltonian (2.1), (2.13), and deduce from this that 
their spectra coincide. The degeneracies are also explained in terms of the 
degeneracies of the N = 2 system (see also Sec. 6). 

Since the Hamiltonian (2.1) commutes with the operator Z (2.18), we can 
choose its eigenstates among the eigenstates of Z, i.e., with definite value for 
the number of states "0" (and hence for the number of states "1"). The 
Hamiltonian (2.15) also commutes with Z, and moreover it commutes also 
with all the J2i E™ m for to = 1, ■ ■ ■ , N — 1. All the charges "to" are then 
conserved. It is actually easy to check that the eigenstates of (2.1) are also 
eigenstates of (2.15). They define a set of 2 L eigenstates of (2.15) that involve 
only the states "0" and "1". So we already know that the spectrum of (2.1) 
is included in that of (2.15). 

We will now introduce a set of operators that commute with the Hamil- 
tonian (2.15), and then use them to construct all the eigenstates of (2.15) by 
action on the previously defined eigenstates that involve only the states " 0" 
and"l". 

Let us define 

J ( ™ } ; n,TO = l,...,7V-l; r = l,---,L; s = l,---,r (3.1) 
as the operator with the following properties: 

a) J™ s) Z = ZJ^ s) = (L — r)Jf r m s) , i. e., J ( n ™ } vanishes on states with a 
number of "0" different from L — r. 

b) acts as E nrn on the s^ 1 site on which the state is different from 
"0". (Note that s < r and that we count from the left to the right.) 

For instance, 

Jf 4 2 3) |201201) = |201301) 

J^) |201201) = if r ^4. (3.2) 

We then have 

EEE^M = !- J (o)^ (3-3) 

r=l s=l m=l 
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where J(o) is the projector on the state z — L, i. e., with only "0". By 
straightforward calculations one can show that the Hamiltonian (2.15) com- 
mutes with all the J£™ (r = 1, • • • , L ; s = 1, • • • , r ; n, m = 1, • • • , N — 1) 
and J( ). This property is only valid for the open chain. 

Consider now an eigenstate ip of the Hamiltonia (2.15), involving only 
the states "0" and "1", and let r be the eigenvalue of L — Z on it. Then the 
wave functions defined by 

n j;i v> (3.4) 

s=l 

are eigenstates of the Hamiltonian (2.15) with the same energy as ip. Here 
the n s take the set of values 1, • • • , N — 1 for s = 1, • • • , r. The set of these 
wave functions with (ni, • • • , n r ) G {1, • • • , N — l} r is linearly independent. 

So each ifj eigenstate of (2.1) leads to (N — l) r eigenstates of (2.15). The 
number of ^'s such that (L — Z)ip = rip being (^j , we get a total of 

(N — l) r = N L (3.5) 

eigenstates. Thus all the eigenstates are given by Eq. (3.4) and the spectrum 
of (2.15) is identical with that of (2.1). This construction also provides the 
degeneracies of eigenvalues of (2.15) in terms of those of (2.1): each state of 
a multiplet of (2.1) provides (N — l) r states of the corresponding multiplet 
of (2.15). We conclude this section with the following observations: 

a) This construction does not work for periodic chains. (The do not 
commute with the Hamiltonian in this case.) 

b) As explained in the following sections in the example of iV = 3 the 
Hamiltonian (2.15) is in fact invariant under U q (sl(N — 1)) for any q. 
To each site is associated the sum of the trivial ("0" state ) and the 
fundamental ("1" to "N — 1" states) representations. The previous 
construction can be seen as a consequence of the fact that the iterated 
tensor products of the fundamental representation of U g (sl(N — 1)) are 
irreducible [13] (See Appendix C). 

c) The A^-state generalization of the Hamiltonian with the constraints 
(2.5) has the diagonal form given by Eq. (2.6) where now n^{k = 
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1, • • • , L — 1) takes the values one (one time) and zero (N — 1 times). 
This is a generalization of the fermionic occupation number. 

4 Invariance of the Hamiltonian 
under U q (sl{2)) for all q 

In the present and in the subsequent sections, we will discuss the relation- 
ship between the spectra of the 2-state Hamiltonian (2.1), (2.13) and the 
3-state Hamiltonian (2.15), (2.16) from a different point of view, namely, the 
[/g(s/(2))-invariance of the latter. As indicated, we choose N = 3, although 
our discussion holds in the general case, too. As we first noticed by numer- 
ical diagonalizations, the spectrum of the 3-state Hamiltonian H is highly 
degenerate. This suggests that the symmetry of H should be much higher 
than the obvious SU(2) ® U(l) invariance mentioned in the introduction. 
Visibly, this invariance is generated at one site by the matrices e^, e z , e° 
specified in (1.7). Thus the 3-dimensional state space at one site carries the 
direct sum of a singlet and a doublet representation of SU (2) (corresponding 
to the indices and 1,2, respectively). 

To gain some hints towards the symmetry of H we sat about to construct 
a g-analogue of this Hamiltonian. This is easily done, once one has noticed 
that the doublets (t^ and (t ;± ) defined by 

t + = E 10 , r = E 20 (4.1) 

t' + = E 02 , t'~ = -E m (4.2) 
are tensor operators of SU{2) and that, up to the normalization, 

(t ® t') = -E 10 ® E 01 - E 20 ® E m (4.3) 

(t' ® t) = E 02 ® E 20 + E 01 ® E 10 (4.4) 

are the scalar parts of the tensor product of t and t' and of t' and t, respec- 
tively. Notice that in Eq. (1.5) we put t~ = t + , p~ = t~ , p + = t' + and 

To g-deform this structure we follow Ref. [14]. We recall some of the def- 
initions. Let q be a non-zero complex number with q 2 ^ 1. Then U q (sl(2)) is 
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the universal algebra (associative, with unit element) with generators e, /, k ±l 
and relations 

kk- 1 = k~ l k = 1 (4.5) 
kek~ l = qe , kfk~ l = g" 1 f (4.6) 
P - P 2 

ef~fe= (4.7) 
q-q 1 

(e and / correspond to the generators S' ± described in Sec. 1). It is converted 
into a Hopf algebra by introducing the coproduct A through 

A(e) = e <g> k~ l + k®e 

A(/) = f®k~ l + k®f (4.8) 
A(P^) = k ±l ® fc ±1 , 



the co-unit e through 



5(e) = £ (/) = 

e(k*) = 1, (4.9) 



and the ant ip ode 5 through 



5(e) = -g^e 

^(/) = -?/ (4-10) 
S(P^) = fc=Fi. 



It is well known that the spin \ representation of U q (sl(2)) is undeformed 
(and, of course, so is the spin representation). Thus the relevant represen- 
tation of U q (sl(2)) at one site is undeformed as well and is given by 

e _> E 12 , / — > £ 21 , fc -> £ 00 + g 1 / 2 ^ 11 + g- 1 / 2 ^ 22 . (4.11) 

Consequently, the deformation of the representation of U q {sl{2)) in the state 
space of H is introduced by the non-trivial coproduct. 

The g-deformed tensor operators (t^) and (t'^) are easily found to be 

t+ = t + = E w , t~ = t~ = E 20 (4.12) 
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t'+ = t' + = E 02 , t'- = qt'- = -qE m . (4.13) 

Using the techniques of Ref. [14], we can then calculate the g-deformed 
versions of Eqs. (4.3), (4.4). The surprising result is that with an adequate 
normalization 

(t q ®t' q )o = (W)o (4.14) 

(t q ®t g ) = (t'®t) (4.15) 

(the underlining of <8> is to indicate that the tensor product is constructed 
according to the rules for tensor operators over quasitriangular Hopf alge- 
bras). Thus the scalar operators in Eqs. (4.3), (4.4) are undeformed. This 
is not completely trivial since, again with a suitable normalization, we find 

(t q ®t q ) = E 10 <g> E 20 — qE 20 (g) E 10 (4.16) 

(t' q ®t' q ) = E 01 (g) E 02 — qE 02 <g> E 01 ; (4.17) 

thus, these scalar operators are deformed. 

Since the remaining terms of H are obviously ?7 g (s/(2))-invariant, we con- 
clude that the undeformed Hamiltonian H is f/ g (s/(2))-invariant, for all q. 
In view of the well-known homomorphism of the g-deformed affine algebra 
U q {sl(2)) onto U q (sl{2)) (see the subsequent section), it comes to the same 
to say that H is even f/ g (s/(2))-invariant, for all q. 

5 Construction of ?7 g (5/(2))-invariant 
Hamiltonians 

The upshot of the foregoing section was that the Hamiltonian H discussed 
there is t/ g (s/(2))-invariant. In the present section we are going to rederive 
this result by means of a more systematic approach. We shall start from a 
class of Hamiltonians generalizing H and derive the conditions under which 
they are f/ g /(s/(2))-invariant for a fixed value q' of q. This will lead us back 
to the Hamiltonian H . 

To begin with, let us recall the definition of the Hopf algebra U q (sl(2)). 
With q chosen as above, this is the universal algebra (associative, with unit 
element) generated by e i: fi, kf 1 ; % E {1,2}, with relations G {1,2}): 

Jv^Jv^ — — Jv^ Jvi — — X , kjkj fvjfvi ^5.1) 
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kiejki 1 = q^ej , hfjK 1 = Q fj 



(5.2) 



k 2 - k~ 2 

e ifi ~ fi e i = $ij ~ _ \ (5-3) 



e\ej - [3],e?ej-ei + [^e^ 2 - e^-ef = 0, i^j (5.4) 

./;'./} - [3],/?/;/* + [3] q fifjff - fjff = 0,i^3, (5-5) 

where 

(«,') = (_V 2 2 ) (5.6) 

is the Cartan matrix of the affine Lie algebra sl(2) and where, for any integer 
n, we define 

M, = ^f . (5.7) 
The Hopf algebra structure is defined by introducing the coproduct A through 

A(ej) = e { ® k^ 1 + k { <g> 

A(/i) = /i ® fcf 1 + h<S> fi (5.8) 
A(fcf) = fcf 1 ® fcf 1 , 



the co-unit e through 



e(e t ) = eifi) = 

e(k? 1 ) = l , (5.9) 



and the antipode 5 through 

= (5.10) 
^(fcf 1 ) = fcf- 
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For our purposes it is important to recall that for any non-zero complex 
number c there exists a unique algebra homomorphism g of U q (sl(2)) onto 
U q {sl{2)) such that 

g(ei) = e , g(e ) = cf 

9(fi)=f , g(fo)=c- 1 e (5.11) 

g(h) = k , g(k ) = k' 1 . 

Consequently, if u is any representation of U q (sl(2)) in a vector space W, then 
uj o g is a representation of U q (sl(2)) in that vector space. The coproduct of 
U q (sl{2)) can then be used to construct the L tensorial power of uog. We 
stress that this representation is different from that obtained by composing 
the tensorial power of cu with g [for g is not a coalgebra homomorphism, 
i.e., in general we have Ag(a) ^ (g <g> g)A(a) for a G U q {sl{2)) and hence, for 
example, (u; ® u;)Ag(a) 7^ (a;^ <8> a;<?)A(a)]. Visibly, the representation o g 
depends on c. In the following we shall choose c = 1. 

After these preliminaries we proceed to our main topic. Let A be any 
finite index set not containing the numbers 1 and 2, and let n be the number 
of its elements. Later in this section we shall choose A = {0}. In the 
following, the Greek indices a, (3, 7, 5 run through A U {1, 2}, and the Latin 
indices a, b, c, d through A. 

We consider the 2-sites Hamiltonian 

H 12 = Y,c a ^sE a P®E^ s , (5.12) 

where the c a/ 3 y s are arbitrary complex numbers, the E a/3 are the basic 
(n + 2) x (n + 2) matrices operating on W = C ra+2 (C denotes the field 
of complex numbers), and where if 12 is an operator in W ® W. We choose 
a non-zero complex number q' with g /2 7^ 1 (later on we shall assume that q' 
is not a root of unity). According to the foregoing discussion there exists a 
unique representation p of U q i(sl(2)) in such that 

p( ei ) = p(/ ) = £ 12 

= P(eo) = E 21 (5.13) 
p(fci) = p{K v ) = K , 

with 

K = J2E aa + q' 1/2 E 11 + q' ~V 2 E 22 . (5.14) 
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Using the coproduct of U q r(sl(2)) we obtain the tensorial square (p ® p) o A 
of p, i. e., the following representation of U q /(sl(2)) in W <S> W: 

{p®p)A{e 1 ) = E l2 ®K~ l + K®E 12 

(p<g>p)A(/ ) = E 12 ®K + K~ l ®E 12 

(p®p)A(f 1 ) = E 2l ®K- x + K®E 21 

(p®p)A(e ) = -E 21 ® If + K" 1 ® _E 21 (5.15) 

(p®p)A(A;i) = K ® K 

{p®p)A{k Q ) = K- X ®K- X . 

One can show that if 12 commutes with (p <g> p)A(ej) and (p <g> p)A(fj); j G 
{1,2}, if and only if it takes the following special form: 

#12 = J2 C abcdE ab ®E Cd 
abed 

+ J2 c abiiE ab ^(E u + E 22 ) 

ab 

+ c ucd(E u + E 22 ) ® E cd 

cd 

+ E Caiid(E al <g> E ld + £ a2 <g> £ M ) 

ad 
be 

+ ci m (£ u + £> (£ n + £ 22 ) . (5.16) 

We stress that the parameters appearing in this Hamiltonian are arbitrary 
complex numbers. Thus, no (/-dependence has been introduced by our re- 
quirement for £/q'(sZ (2) )-invariance. Consequently, the Hamiltonian (5.16) is 
C/ 9 (sZ(2))-invariant for all q. 

It is easy to see that H 12 has the following more general symmetries. 
Define 

D = q J2E aa + qiE u + q 2 E 22 (5.17) 

a 

D' = q J2E aa + q[E u + q' 2 E 22 , (5.18) 

a 

where ?0j 9ij 9ij 92, ?2 are arbitrary complex numbers. Then H\ 2 commutes 
with 

E 12 <g> D + D' <g> E 12 , E 21 ® D + D' ® E 21 , D&D, (5.19) 
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and, of course, also with D' <g> D' . Obviously, this implies that the corre- 
sponding L-sites Hamiltonian 

L-l 

H = ]T Hi, i+ i (5.20) 
i=i 

commutes with the corresponding symmetry operators in W® L , i.e., with 

i 

^D'^^B 21 ®]) 8 ^, (5.21) 

i 



The relevance will be illustrated in the following example. We shall restrict 
our attention to the case A = {0}; furthermore, we shall assume that c 0110 
and ciooi are different from 0. In that case, we may even assume that 

c ono — c iooi — — 1) (5.22) 

for up to the overall normalization of H this can be achieved by a similarity 
transformation in the state space W® L (recall that W = C n+2 = C 3 ), which 
preserves the special form of H. With (5.22), H takes the form (2.15), (2.16), 
with 

V = Clin — Cnoo — Cooil + Coooo 

w = -(cnoo + coon J - coooo 
a = 2( Cl100 — c oon) 

b = coooo • (5.23) 

Thus, choosing A = {0} and assuming that (5.22) is satisfied, the U q >(sl(2))- 
invariant Hamiltonians of the form (5.12) are exactly the Hamiltonians in- 
troduced in (2.15), (2.16). 
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It is easy to see what the symmetry transformations (5.21) imply for the 
spectral decomposition of H. Let vo, i>i, i>2 denote the canonical basis vectors 
of W = C 3 . For any subset J C {1, 2, • • • , L}, let Vj denote the subspace of 

V = W® L , (5.24) 

which is spanned by the tensors v ai <S> ■ • ■ <S>v aL , with a, G {1, 2} if i G J, and 
a, = otherwise. Obviously, V is the direct sum of the subspaces Vj and 

dimVj = 2 r with r = p. (5.25) 

Furthermore, the Vj are invariant under the transformations (5.21), and it 
can be shown that these transformations act irreducibly on Vj. 

th 

The latter result can also be derived as follows. Visibly, Vj carries the r 
tensorial power of the elementary spin \ representation of U q >(sl{2)). But if 
q' is not a root of unity, it is known [13] that this representation is irreducible. 

The upshot is that V = W® L decomposes into irreducible 2 r -dimensional 
subspaces, with r e {0, 1, ■■-,£,}, and that exactly (^j of these subspaces 
have that dimension. This implies that the eigenspaces of H also decompose 
into irreducible subspaces of these dimensions. 

To be more precise, let r e {0, 1, • • • , L} and let 

V r = Vj. (5.26) 

Thus V r is the eigenspace of the operator Z (defined in Eq. (2.18)) corre- 
sponding to the eigenvalue L — r. We stress that for r fixed, the irreducible 
representations acting in the Vj with %J — r are equivalent. 
Furthermore, for a G {1, 2} we define 

V a = (Cv + Cv a f L , (5.27) 

i.e., V a is spanned by the tensors of the form x\ <S> • • • <S>xl, with Xi G {i>o, v a } 
for all i. Finally, let 

V r a = V a n V r . (5.28) 

Obviously, we have 

L 

V a = V r a . (5.29) 

r=0 
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Furthermore, it is easy to see that H maps the spaces V a and V r into them- 
selves, hence also the spaces Vf- Thus, according to the foregoing discussion, 
we only have to construct the spectral decomposition of H in the subspaces 
V® for one value of a, say a — 1. The action of the symmetry transforma- 
tions will then yield the corresponding decomposition of V r , hence of V. In 
particular, if the multiplicity of an eigenvalue of H in is equal to to, then 
the multiplicity of that eigenvalue of H in V r is equal to to • T . 

Summarizing the preceding discussion we see that by means of the symme- 
try transformations the spectral decomposition of H in V has been reduced 
to that of H in one of the subspaces V a , i.e., to that of the 2-state Hamilto- 
nian (2.1), (2.13). This is in perfect agreement with the results of Sec. 3. In 
particular, the existence of the symmetries described in Sec. 3 follows 
directly from the results of the present section. This will be explained in 
Appendix C. 

6 Multiplets of the 3 L chain when 
the underlying Heisenberg chain 
has supplementary symmetries 

Consider the 3 L chain obtained with the construction of Sec. 5 from the 2 L 
chain. We use the notation 

d\ = w — a , d,2 = w + a , d% = v + 2w (6.1) 

for the diagonal parameters. As explained in Appendix B, we have to con- 
sider four cases for the values of d\,d2,d 3 . 

Case la 

di = q~ l , d 2 = q , d 3 = q + q^ 1 . (6.2) 

The underlying 2 L chain then has the L r _ g (s/(1/1)) symmetry, and its states 
build doublets corresponding to adjacent eigenvalues r and r — 1 of L — Z. 
In the 3 L chain, the two vectors of each doublet are highest-weight vectors of 
multiplets for the symmetry U q >(sl(2)), with dimensions 2 r and 2 r_1 . These 
multiplets are gathered in a higher multiplet of dimension 2 r + 2 r ~ 1 = 3-2 r_1 . 
There are (^Zi) multiplets of this dimension. In the fermion picture, (^Zi) 
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is the number of states with r — 1 fermions (among L — 1 fermions with 
non-zero energy). 

Case lb 

d 1 = q' 1 , d 2 = q , d 3 = . (6.3) 

The underlying 2 L chain has U q (sl(2)) symmetry, and the states are gath- 
ered into multiplets corresponding to the decomposition of (Spini) (Xli with 
the fusion rules of U q (sl{2)). We shall consider only the case where q is not a 
root of unity. [The case where it is would easily follow from the decomposi- 
tion of (Spin|) 0L into irreducible and indecomposable representations.] The 
multiplets of the 2 L chain collect states with different eigenvalues of L — Z 
as follows: 

{L,L-1,...,0} , [L + l] , (1) 
{L-l,L-2,.-.,l} , [L-l] , (L-l) 

{r,r-l,--.,L-r} , [2r - L + 1] , (( L L _ r ) - ( L ^_ 1 )) 

with 2r > L. We have put in braces the eigenvalues of L — Z. The brackets 
indicate the dimension of this L r q (s/(2))-multiplet, whereas the parentheses 
indicates the number of multiplets of this type. 

For the 3 L chain, each state of these L^(s/(2))-multiplets becomes a 
highest-weight vector of a L^/(sZ(2))-multiplet. We then have larger mul- 
tiplets of dimensions 

2 l + 2 l-i H h 2° = 2 L+1 - 1 , (1) 

2 l-i + 2 l- 2 + ... + 2 i = 2 L -2 , (L-l) 

2 r + 2 r - 1 + ... + 2 L - r = 2 r+1 -2 L ~ r , (( L L _ r ) - (l^)) 
where again, in the parenthesis, we indicate the numbers of these multiplets. 

Case 2 

d 1 d 2 = l , (di-d 3 )(d 2 -d 3 )^l. (6.4) 

As shown in Appendix B, the supplementary symmetry builds a unique 2- 
dimensional multiplet of the 2 L chain with the state v® L [for the notation, 
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see the lines preceding Eq. (5.24)] and a state with the eigenvalue of Z equal 
to L — L In the 3 L chain, these two states become highest-weight vectors 
of L^/(s/(2))-multiplets of dimension 1 and 2, respectively. So we finally get 
a unique multiplet of dimension three, while the other [/ g /(s/(2))-multiplets 
are not combined. 

Case 3 

did 2 7^ 1 , (di - d 3 ) (d 2 - d 3 ) — 1 . (6.5) 

The two states that build a doublet in the 2 L chain have now Z-eigenvalues 
and I. In the 3 L chain problem, these states are highest-weight vectors of 
f/ g /(s/(2))-multiplets of dimensions 2 L and 2 L ~ 1 . We then have a multiplet 
of dimension 3 • 2 L_1 . All the other states remain in multiplets of dimension 
T (0 < r < L). 

7 Conclusions 

We have started our paper with a physically interesting one-dimensional 
system of fermions with an interaction given by three coupling constants 
Y v Y t and ^. The phase diagram of this system is unknown. One rigorous 
result of our investigation is that in the strong repulsion limit, the phase 
diagram coincides with the one of the Heisenberg chain in a magnetic field 
(see Fig. 1). The Eq. (2.14) of Sec. 2 describes this mapping. We also give 
an interpretation of the phase diagram of Fig. 1, which is usually understood 
in the language of magnetic systems, in terms of hard-fermion physics (see 
Appendix A). Our work has gone beyond this investigation. On the one 
hand in Sec. 2 we have obtained TV-state models whose spectra coincide with 
those of the Heisenberg model (the N = 3 case being relevant to the fermionic 
problem). The eigenfunctions of these Hamiltonians can be obtained from 
those of the N = 2 case using the construction given in Sec. 3 (see also 
Appendix C). On the other hand we have noticed (see Sec. 4) that for 
N = 3 the multiplet structure of the chains is related to finite-dimensional 
representations of the affine algebra U q (sl(2)), the chains being invariant 
under this algebra for any q. In Sec. 5 we give the conditions under which 
an (n + 2)-state Hamiltonian is invariant under U q {sl{2)) for any q and find a 
new set of Hamiltonians whose spectra are no longer given by the Heisenberg 
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model. Finally, we would like to notice that for the Heisenberg model we have 
found a new symmetry valid for finite chains when the parameters verify the 
relations (2.11). This symmetry is discussed in Appendix B. In particular, 
this symmetry applies to the Pokrovsky-Talapov line. 
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Appendix A From 4-state fermionic to 

3-state magnetic spin chains 

In this appendix we obtain the Hamiltonian (1.3)— (1.5) as the U — > oo limit of 
the fermionic Hamiltonian (1.1), (1.2). Since in this limit double occupancy 
of fermions is forbidden (it costs an infinite energy to put fermions at the 
same site) we can start our derivation by constructing a Fock space where, 
in a given lattice point, we may have no fermions ("0"), a fermion with up 
spin ("T") or down spin ("J."). Due to the fermionic nature of the operators 
appearing in (1.3)— (1.5) we choose an order and construct directly a repre- 
sentation in the 3 L -dimensional Fock space. The chosen order is illustrated 
in the following two examples 



T,T,l,o,t,...,T) = c+ T c+ T c3+ iC + T ...c+ T |o,o,o,...,o) 
|o,t,|,o,t,...) = c+ T c^c+ T ... |0,0,0,...,0), 



where 



CV|0,0,0,...) = 0, i = l,2,...,L; a =U • 
With this convention we can represent the Hamiltonian (1.1), (1.2) as 



(A.l) 
(A.2) 



L-l 

H' = £ li 

i=i 



l»-i< 



d 1 \-l 
d 1 \ 

-1 i do 



-1 



-1: 



da 



!l (A.3) 
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r> L-i 
H[ = - E li 

1 *tk 



-1: 




-1 



-2 



(A.4) 



where d = 0,^ = — ^,rf2 = ^- + ^,^3 = 7 + 7- and 1, is a 3 x 3 unit 
matrix attached at site i. The elements not shown explicitly in the above 
matrices have zero value. 

This Hamiltonian can promptly be rewritten in terms of the matrices in- 
troduced in (1.7) 



l-\ r 

H 'o = El - (ft ft+i + QiQT+i + T ? T r+i + 
i=i 

+ y«n + | {4 + + ^ ( e ° - + 6} (A.5) 

H[ = -%j:el. (A.6) 
zr i=i 

The reader has already noticed that our proof of the connection between the 
Hamiltonians (1.3)— (1.5) and (1.1) and (1.2) applies to open chains only. In 
the case of a periodic chain one has to split the sectors of the fermionic Fock 
space according to the parity of the fermionic number. Sectors with an odd 
(even) number of fermions will correspond to the Hamiltonian (1.3)— (1.5) 
with periodic ( ant iper iodic) boundary condition. 
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Appendix B New symmetries of 

the Heisenberg model 

Let us consider the 2 L state Hamiltonian H = J2j=i Vj with Vj given by Eq. 
(2.13); Vj can also be written 



Vj = b + 



\ 



di :-l 
-li d 2 



V 



with d\ = w — a , d 2 = w + a , rf 3 = f + 2u>. 
We look for symmetries of H of the form 



(B.l) 



3=1 



B.2) 



where -D and D' are diagonal matrices and D®^ -1 ) = D D . . , 
(j — 1 times). 

If if commutes with e given by (B.2), it will also commute with /, its 
transpose, given by 



(B.3) 



The existence of these symmetries depends on the values of di,d 2 ,ds. Let us 
write 



D 



(B.4) 



Then the conditions on l i: l' i ,d 1: d 2: d 3 for [H, e] = are 

cMd -l = , 1' - d 2 l = 
{d 1 - d 3 )h -l[ = , h- (d 2 - d 3 )l[ = . 



(B.5) 
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The solvability of these equations depends on the values of di, d 2 , d 3 , and we 
have to consider four cases: 

Case 1 

d 1 d 2 = l , {d 1 -d z ){d 2 -d z ) = \ (B.6) 

i.e., either 

di = q' 1 , d 2 = q , d 3 = q + q' 1 (B.7a) 

or 

di = g" 1 , d 2 = q , d 3 = . (B.7b) 

These cases correspond to the conditions (2.5) and (2.9), respectively, and 
the symmetries generated by e and / can be identified with U- q (sl(l/l)) and 
U q (sl(2)), respectively. 

Case 2 

d l d 2 = l , (d 1 -d 3 )(d 2 -d 3 ) ± 1 . (B.8) 

Then 

^(o°o). H^o)' (B ' 9) 

Dropping an overall factor Zq -1 , we have 

e = d{-\E™)^-V ® E 01 <g) (E 00 )®^-^ (B.10) 
i=i 

/ = jr di-\E 00 )®^ E io (£00)®(i-i) (B n) 

We then see that all the states are cancelled by /, except the state with 
charges "0" only. In fact, e and / connect the states |0 . . . 0) and /|0 . . . 0) 
in a doublet, and all the remaining states are singlets. 

Two levels are then degenerate, the others being non-degenerate. This 
degeneracy of two levels only has also been checked numerically. 

Case 3 

d x d 2 ^\ , (d 1 - d 3 )(d 2 - d 3 ) = 1 (B.12) 
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The situation is similar to the previous case and we can thus write 

e = jr ( dl - d,y- 1 (E^Y^ ® E m ® (E^Y [L - j) (B.13) 
j'=i 

/ = £ (di - dsY' 1 (E^Y U - 1] ® E 10 ® . (B.14) 

Now e vanishes on all the states but |1 . . . 1), and {|1 . . . 1), e|l . . . 1)} build 
a doublet. All the other states are singlets for e, /. Again, only two levels 
coincide. 

Case 4 

d t d 2 ± 1 , (d 1 - d 3 )(d 2 - d 3 ) ± 1 . (B.15) 
Then e = / = 0; we get no further symmetry. 

The symmetries that exist for the 2 L state Hamiltonian in the cases (1), (2), 
and (3) are combined with the other symmetries [U q >(sl(N)) in the (iV + 1) L 
state case]. In particular, when two states build a doublet of the 2 L state 
Hamiltonian, the representations of U q i(sl(N)) associated with these states 
m the (N + l) L -state case are also combined into higher multiplets. In the 
case (2), we get a 1 + iV-dimensional multiplet, together with iV r -dimensional 
multiplets. In the case (3), we get a N L + N L ~ ^dimensional multiplet, to- 
gether with iV r -dimensional multiplets (see Sec. 6). 

Appendix C The symmetries J^) 

reconsidered 

In the present appendix we want to show that the existence of the symmetries 
J™!s) described in Sec. 3 follows from the results of Sec. 5 by using the 
standard representation theory of associative algebras. In doing so we use 
the notation introduced in Sec. 5. 

Let p L be the tensorial power of the representation p of U q (sl(2)) as 
specified in Eq. (5.13) (we simplify the notation and write q instead of q'). 
Thus p L acts in the state space V = W® L of the L-sites chain, and we know 
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that V decomposes into the direct sum of the invariant irreducible subspaces 
Vj, J C {1,2, ... , L}. Obviously, if J, J' C {1,2,..., L} and if %J and fjj' 
are different, the representations given by Vj and Vj> are inequivalent. On 
the other hand, if JjJ = %J' = r, these representations are equivalent. More 
precisely, if we forget about the tensorial factors v (which can be done by 

use of canonical ?7 g (s/(2))-module isomorphisms) they simply coincide with 
th 

the r tensorial power of the elementary 2-dimensional representation of 
U q (s7(2)), acting in (Cui + Cv 2 )® r . 

Now the representation theory of associative algebras tells us the follow- 
ing. Choose r e {0, 1, . . . , L} and let A be an arbitrary linear operator in 
(Ci>i + Cv2)® r ■ Then there exists an element a E U q (sl(2)) such that p L (a) 
acts on Vj as A if (j J = r, but as if jj J ^ r. 

Of course, all these operators p L (a) commute with the Hamiltonian H; 
moreover, it is easy to construct the operator p L (a) explicitly. Without loss 
of generality, we may assume that 



where Ax, . . . , A r are arbitrary 2x2 matrices, with matrix indices in {1,2}. 
Let A s be the 3x3 matrix obtained from A s by adding a zeroth row and 
column consisting of zeros only. Define the operator J r (Ai, . . . , A r ), acting 
in V, as follows: 

J r {M, ■ ■ ■ , A) = E °° ® • • • ® Ai ® ■ ■ ■ ® M ® ■ ■ ■ ® Ar ® ■ ■ ■ <S> E 00 , (C.2) 

where the sum is extended over all tensor products of the matrices A ± , . . . , A r 
and L — r matrices E 00 , with the Ai, . . . , A r appearing in their given order. 
Then J r (A\, . . . , A r ) is the operator acting on Vj as A\ <g> . . . <g) A r if JJJ = r 
and as if ft J ^ r. Thus we know that the operator J r (Ai, . . . , A r ) commutes 
with the Hamiltonian H. 

We note that if r' G {0, 1, . . . , L} and if A' 1: . . . , A' r , is a second family of 
2x2 matrices, we have 



Moreover, our results show that the operator algebra p L (U q (sl(2))) coincides 
with the linear span of the set of all operators J r (Ai, . . . , A r ). In particular, 



A = A x ® ■ ■ ■ ® A, 



(C.l) 



J r (Ax, . . . , A r ) JyiyA'-y, . . . , A' r ,) — J r yAxA' 1: . . . , A r A' r 

= 



) if r — r' 
if r^r'. 



(C.3) 
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the operators p L (ei), p L (fi), p L (kf r );i E {1,2}, can be written in this form. 
For example, by definition of p L , we have 

P L (ei) =EE Mk, E 12 , k~\ k- 1 ) , (C.4) 

r=l s=l 

where -E 12 is in position s and k denotes the 2x2 matrix 

k = q 1/2 E u + q~ l/2 E 22 . (C.5) 

On the other hand, it is easy to see that 

JJ5 = J r (l,...,£™\...,l), (C.6) 

where, now, E mn is in position s and 1 denotes the 2x2 unit matrix. In 
view of Eq. (C.3), this implies that the operators and J( ) generate the 

associative operator algebra p L (U q (sl(2))). 
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